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Density-wave patterns in (quasi-) discrete media with local interactions are known to be unstable. 
C\ . We demonstrate that stable double- and triple- period patterns (DPPs and TPPs), with respect to 

D ' the period of the underlying lattice, exist in media with nonlocal nonlinearity. This is shown 

[JL^ ' in detail for dipolar Bose-Einstein condensates (BECs), loaded into a deep one-dimensional (ID) 

optical lattice (OL), by means of analytical and numerical methods in the tight-binding limit. 
The patterns featuring multiple periodicities are generated by the modulational instability of the 
continuous-wave (CW) state, whose period is identical to that of the OL. The DPP and TPP 
C/3 , emerge via phase transitions of the second and first kind, respectively. The emerging patterns may 

be stable provided that the dipole-dipole (DD) interactions are repulsive and sufficiently strong, in 
comparison with the local repulsive nonlinearity. Within the set of the considered states, the TPPs 
realize a minimum of the free energy. Accordingly, a vast stability region for the TPPs is found in 
C , the parameter space, while the DPP stability region is relatively narrow. The same mechanism may 

create stable density-wave patterns in other physical media featuring nonlocal interactions, such as 
arrayed optical waveguides with thermal nonlinearity. 
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PACS numbers: 03.75.Lm; 05.45.Yv 



Introduction. The formation of density-wave patterns in periodically structured media is a fundamental physical 
effect, celebrated manifestations of which are the Peierls instability of the electron gas in one- dimensional (ID) metals, 
leading to the emergence of charge-density waves P, Q , and spin-density waves [2] . Recently, many complex settings 
^ ' known in solid-state physics have been reproduced {simulated) in an essentially simpler form in quantum gases — 
'—'I often, these are atomic Bose-Einstein condensates (BECs) — loaded into optical lattices (OLs) Q. Thus far, this 
^ . possibility was not demonstrated for density-wavepatterns, the reason being that such patterns formed in EEC with 
^ ' contact interactions between atoms are unstable [J]. The objective of this work is to demonstrate that stable density 
ly-j , waves are possible in BEG with long-range dipole-dipole (DD) inter-atomic interactions. The same mechanism should 
make it possible to create stable density (intensity) waves in other periodically structured physical media featuring a 
T-H nonlocal nonlinear response, which is inherently present in the transport mechanisms of heat and charge carriers 
■ 6|, electrostatic interactions in liquid crystals 7], photon self-attraction and many-body interactions in plasmas 
psj i m ^'^d BEG (Iflj. Further, an effective gravitation-like attraction between atoms can be optically induced in BEG 
i by means of resonant illumination It was demonstrated theoretically and experimentally [11] that the 

' implementation of patterns supported by the nonlocal nonlinearity is especially relevant in thermal optical media, 
, where the periodic structure can be readily built in the form of arrayed waveguiding systems. 
J> ' BEG of dipolar atoms or molecules loaded into OLs have been in the focus of research work for the past decade, 
; ""j ', startin g fro m the theoretical analysis of basic properties of such condensates [l^ , followed by the analysis of quantum 
rS phases [15l.[l6j. various textures [13, and suppression of the quantum collapse (fall onto an attractive center) [l^ 
in them. Structured ground states and supersolidity of dipolar gases [3, [13, [HI , as well as the rotonic dispersion 
" " " relation [l^], were discussed theoretically but still await experimental confirmation. Experimentally, first effects of 
dipolar interactions in quantum gases were observed in the BEG of ^^Gr atoms [1^. A condensate of ^^^Dy atoms 
was recently created too [IJ, which provides an even stronger DD interaction than ^^Gr. In addition to that, the 



creation of dipolar BEG is expected in erbium [25| and in gases of molecules carrying electric dipole moments |26| . 
A powerful tool which allows one to control dynamical effects in the dipolar BEG, such as the onset of collapse, is 
the use of the Fcshbach resonance for tuning the strength of contact interactions between atoms, which compete with 
their long-range DD interactions [27j . The results obtained in this field have been summarized in recent review [28| . 

Another ubiquitous tool used in experiments with dipolar BEG is provided by OLs. Recently, dipolar effects in the 
chromium condensate trapped in OLs were reported in Refs. [29|-|3]|. In particular, strong DD repulsion can stabilize 
the condensate with attractive contact interactions, trapped in the OL [3i|. It is well known that the OL facilitates 
the creation of gap solitons [32], or bound pairs of repulsively interacting atoms, at the microscopic level |33l] . Gap 
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solitons in the ID model of dipolar condensates trapped in the OL potentials were studied in Ref. '34']. 

While the state of the trapped condensate may obviously feature the same periodicity as the underlying lattice (we 
call it the continuous- wave state, CW), we here report stable density waves in the form of double- and triple-period 
patterns (DPPs and TPPs) in the dipolar BEC loaded into a deep OL potential. The TPP realizes a relatively deep 
energy minimum in the space of the considered states, and its stability area is much broader than that of the DPP, 
therefore the TPP is a candidate to the role of the ground state. To the best of our knowledge, this is the first 
demonstration of stable density waves of the TPP type in OL-trapped BECs, as well as in dynamical lattices with 
long-range interactions of any physical nature. Featuring long-range off-diagonal and diagonal order different from 
that of the underlying lattice potential, both the DPPs and TPPs may be viewed as a metastable supersolid state 
[20I [21I [sslo f the dipolar condensate. 

In Ref. Q it was demonstrated that the period-doubling instability of the CW in the BEC with local interactions, 
trapped in the deep OL, gives rise to DPPs, but, as mentioned above, the emerging pattern is never stable. As for 
the TPPs, an incentive for their consideration is provided by the recent analysis of the trapped dipolar BEC, which 
concluded that a three- well system is the minimal setting which is necessary to let the nonlocal character of the DD 
interactions manifest itself jlSl]. In this work we consider all possible combinations of the repulsive/attractive contact 
(RC/AC) and repulsive/attractive DD (RDD/ADD) interactions, and demonstrate that unstaggered (regular) DPPs 
and TPPs have their stability regions in the RC-(-RDD case, provided that the long-range DD repulsion is sufficiently 
strong in comparison with the local self-repulsion (these findings imply that, in the opposite AC-I-ADD case, the 
staggered [36J counterparts of these patterns are stable in the respective regions too). The analysis reveals that the 
multiple-period patterns emerge when the corresponding modulational instability of the CW sets in, i.e., the respective 
perturbations seed the creation of the patterns. The TPP's stability area is found to be much larger than that of 
the DPP, and the calculation of the free-energy density demonstrates that the TPP realizes a well-pronounced energy 
minimum, in comparison with the CW and DPP states. When the TPPs are unstable, the DD interactions of either 
sign significantly inhibits the instability, suggesting that even unstable TPPs may be observed in the experiment. 

The model. For condensates loaded into a deep OL, the corresponding Gross-Pitaevskii equation can be reduced 
to the discrete nonlinear Schrodinger equation (DNLSE) in the framework of the tight-binding approximation 0, 113] 
(a similar reduction to the quasi-discrete propagation is well known in optics [38]). In turn, the analysis of various 
stationary states supported by the DNLSE may start from the ant i- continuum limit, which corresponds to the chain 
of uncoupled sites fSy]. The same tight-binding approximation, if applied to the dipolar BEC loaded into a deep OL, 
leads to the ID 15, 40] and 2D |^41] DNLSEs which include the long-range DD interaction between lattice sites. We 
here consider the trapped condensate modeled by the scaled ID equation [40| : 



2 



where a — —I and -fl corresponds to the attractive and repulsive contact interaction, respectively, C > is the 
inter-site coupling constant which accounts for the tunneling of atoms between adjacent sites of the lattice, and F > 
or F < is the relative strength of the DD attraction or repulsion, in comparison to the local nonlinearity. In fact, a 
accounts for the effective on-site nonlinearity, which includes both the DD interaction between atoms trapped in a given 
potential well and the contact interaction proper. Time is measured in units of wj^, where ujj_ is the frequency which 
defines the radius of the transverse confinement, a± — ^ fij (matomWx)- C may also be fixed by rescaling, therefore 
basic findings are displayed below for C — 0.8. In this work, we present results for unstaggered patterns. They can 
be made equivalent to their staggered counterparts by transformation /„ — (— l)"e'**'-'*/Tii {'^'i-T} — )• — {c, F}. 

Equation ([T]) does not include the external trapping potential. Under realistic experimental conditions, the actual 
number of OL periods in the trap is > 100. It is straightforward to check that distortion introduced by this potential 
is negligible for the pattern s rep orted below. Undistorted configurations can also be realized in the dipolar condensate 
loaded into a toroidal trap [43 ]. 

For other nonlocal systems, such as optical media with the thermal nonlinearity [lH, the DNLSE differs by the 
form of the kernel accounting for the long-range interaction. Essentially the same analysis as reported below can be 
developed for such models too. 

Stationary solutions to Eq. ([T]) with chemical potential /i are sought for as /„ — J7„exp {—i/it), which leads to the 
equation for real discrete wave function [/„: 



\m 



Un- (2) 



First, a straightforward analysis of the modulational instability of the unstaggered CW (constant) solution, Ul^ = 
fl/ [a — 2C(3)F], where C(3) w 1.202 is the zeta-function, has been performed, using the linearization of Eq. ([TJ for 
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small perturbations. In the most relevant RC+RDD case, the CW is stable in a relatively narrow parameter region, 
see Fig. [TJ 




FIG. 1: The stability diagram for unstaggered CW states with the repulsive sign of both the contact and DD interactions 
((7 = +1, r < 0). In all figures, the lattice coupling constant is C = 0.8. 



Stable double-period and triple-period patterns. The formation of multiple-period patterns is a consequence of the 
modulational instability of the CW states, which represent uniform chains of BEC droplets trapped in the OL. Ac- 
cordingly, periodically modulated patterns may be called condensate-density waves. They were found numerically by 
solving Eqs. ^ with the help of the modified Powell method [13], and also in exact and approximate analytical forms 
(see below). Their stability was explored by computing eigenvalues for small perturbations, using the linearization of 
Eq. ([1]), which was then checked by direct simulations of Eq. ([1]) for perturbed solutions. Close to the anticontinuum 
limit, patterns with all periodicities are long-lived ones, as the instability weakens with the decay of the inter-site 
coupling. This is a generic property of the DNLSE with the contact nonlinearity [36j , which persists in the presence 
of the DD interaction. 

DPPs which are shown in the top row of Fig. [21 with J7„ taking two values, 0i ^ (t>2, can be found as exact solutions 
ofEq. Q: 
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Figure [3ja) demonstrates that both the linear-stability analysis and direct simulations reveal a rather narrow stability 
area of the unstaggered DPPs for the RC-I-RDD sign combination of the interactions. At large ^, the upper boundary 
of the region is determined by the positivity of (t>i(j)2 (i.e., the "unstaggerness" of the pattern) in Eq. ^ with a = -1-1: 
r < — 2/(3C(3)) ~ —0.55. The DPP is created by the supercritical pitchfork bifurcation [43[ (in other words, through 
the phase transition of the second kind) at = 2C [4FC(3) — 2a\ / \2a -\- 3FC(3)], which is found from Eq. ([3]) by 
setting 01 — (1)2. At the bifurcation point, the DPP branch may emerge as a stable one [Fig. [31[b)], or, at larger |F|, 
as an unstable solution, which quickly enters the stability region with the further increase of /x [Fig. |3^c)]. As seen in 
Fig. ^c), in the latter case the CW solution becomes unstable against other perturbations at /i = /ii, but it remains 
stable against the period-doubling perturbations up to /i = /id-. 



DPP 



TPP 

n= 
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FIG. 2: A scheme of the condensate density waves with the double and triple periodicity. (^1,2,3 are amplitudes of the wave 
function at the lattice sites. 



The TPPs are schematically shown in the bottom row of Fig. [2j Equations for the three real amplitudes, which 
are designated in the figure, are derived from Eq. ([2|): 

= -C(02 + 03 - 201) + a<f>l - r0i {b<f>l + a{4>l + 0^)) , (4) 
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FIG. 3: (Color online) The stability diagram of the unstaggered DPPs (a), and TPPs (d) for the RC+RDD interactions 
((T = +1, r < 0). The DPPs (a)-(c) and TPPs (d)-(f) are stable in the gray area, while the dotted curve, fi = /ii(r), separates 
stability and instability regions (on its left and right sides) of CWs. The DPP is generated, via the pitchfork bifurcation, from 
the CW branch by period-doubling perturbations at — /icr(r), which is marked by the solid curve in (a), while the stable TPP 
is generated by the subcritical pitchfork bifurcation at the turning point, fi — fi2 < fia (d)-(f). Amplitudes of the unstaggered 
DPP and TPP are plotted vs. /i in (b),(c), and (e), (f) for F = —2 and F = —5, respectively [these values of F are marked in (a) 
and (d) by dashed horizontal lines]. Chains of symbols and continuous curves show, severally, numerical results and analytical 
solution for DPP ((S]) and TPP ((5|. The blue line (which seems as a bisector) represents the CW solution. Stable and unstable 
states are shown by solid and dashed lines or symbols, respectively. 
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and two other equations obtained by cyclic transpositions of {1, 2, 3}. Here, a = — [ip{l/3) + ip{2/3)] /54 w 1.157 and 
b = 2C(3)/27 w 0.089, with ?A(z) ee d[ln(r(z))] /dz. These equations give rise to three types of TPPs. Type 1 is 
defined by assuming 4>2 — 4>d. ^ 4>i^ which yields an approximate analytical solution. 



^-2C , , -C a-Vh 



1 = \ TT^ V2 = (P3 = 



(5) 



Type 2 is similar to Type 1, but with comparable absolute values of the amplitudes: > 0, 02 = 03 < 0,and Type 
3 is an exact solution with (pi =0, 02 = —03 = \J {i^ ~ 3C)/ [a — T{a + b)]. The analysis demonstrates that solely 
Type 1 may be stable, therefore only this type is considered below. 

While the unstaggered TPPs exist for all combinations of the repulsive and attractive contact and DD interactions, 
only the RC+RDD combination gives rise to stable states, similar to the DPP. Again, the creation of the TPPs 
is related to properties of the CW background. As seen at Fig. El^d), the CW becomes unstable against generic 
perturbations at /i = //i, but remains stable against the triple-period disturbances up to the bifurcation point, 
M — Mcr, which can be found exactly from Eq. Q: fx„ — — (3C/2) [a — (2a + h)T] / [cr + (a — b)T]. A difference from 
the way the system gives rise to the stable DPP (see above) is that the bifurcation is (weakly) subcritical [i^ in the 
present case (alias it is the phase transition of the first kind), with the stable TPP emerging by a jump at the turning 
point, fj, = IJ.2 < /Xcr, see Figs. [3i;e),(f). 

All the varieties of unstaggered DPPs and TPPs are found to be entirely unstable for other types of the interactions, 
i.e., RC+ADD, AC+ADD, AC+RDD, with the instability growth rate increasing proportionally to the lattice coupling 
constant, C. On the other hand, at fixed C the instability growth rate decreases with the increase of the relative 
strength of the DD interactions (of either sign), which implies that they help to increase the lifetime of unstable 
density waves in the dipolar BEC. Direct simulations confirm the predictions of the linear-stability analysis, as well 
as the inhibition of the instability by the DD interactions. 

A crucial issue is the comparison of energy of coexisting states. Considering spatially unconfined modes, controlled 
by the chemical potential, it is relevant to do this for the mean spatial density (g) of the corresponding free energy, 
G = H — fj,N, where H and N are the Hamiltonian and number of atoms (norm of the wave function). The results, 
obtained analytically for the CW and DPP, and in a numerical form for the TPP, are shown in Fig. U) For a large 
range of parameters the TPP is stable and realizes a relatively deep minimum of the energy, while the DPP may only 
represent a metastable state. Whether the TPP is the ground state, needs to be checked by comparing with more 
complex patterns, if they exist in the system. However even if it is a metastable configuration, it is relevant to explore 
its superfluid properties in further studies. 




FIG. 4: (Color online) The free-energy density of the CW, DPP and TPP states vs. (a) the chemical potential and (b) atomic 
density, Nd, for a = +1 and F = — 5. The results are similar for other negative values of F. 

As concerns parameters relevant to the experimental implementation of the stable condensate-density waves pre- 
dicted here, in chromium the background s-wave scattering length is ^ IOOob. By means of the Feshbach resonance, 
it can be reduced to ~ 2aB, thus making the basic ratio |F| = aDD/os ^ 5, which is sufficient for getting into the 
stability areas displayed in Figs. Elja) and|3l[d). For the dipolar BEC of dysprosium and erbium, the collisional 
properties are not yet known in detail, but similar tuning options are expected. The creation of the DPP may not 
be straightforward, as it does not correspond to the energy minimum, unlike the TPP. To achieve this purpose, a 
periodic modulation of the lattice can be applied, cf. Refs. [3 El]- concerns the observation of the density- wave 
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patterns, in addition to direct imaging, their existence can be confirmed by measuring the momentum distribution in 
the time of flight, after switching off the potential. 

Conclusion. Our analysis predicts stable density waves in OL-trapped dipolar BEG, and, as a matter of fact, in 
other periodically structured media with long-range interactions. In particular, the stable TPPs are predicted for 
the first time, to the best of our knowledge. The stability, which is limited to a relatively narrow region for the 
DPP, and is found in a broad area for the TPP, is provided by the long-range DD interactions, and is impossible 
in the BEG with solely the contact interactions. While the DPP is metastable, the TPP realizes a relatively deep 
energy minimum. These patterns emerge through the phase transitions of the second and first kinds, respectively. 
These condensate-density waves are bosonic counterparts of charge-density waves in ID metals, created by the Peierls 
instability. The results suggest that the long-range DD interactions (whose effective strength in the quasi-lD setting 
can be adjusted by the direction of the polarizing magnetic field) may be used to steer transitions between different 
varieties of Mott-insulator states of the trapped bosons. The same mechanism may generate stable density waves in 
other physical settings — in particular, in optical media with nonlocal nonlinearities. 

The findings reported here suggest a number of potentially interesting developments. A straightforward question is 
whether condensate density waves with periods larger than three may be stable. While this issue should be addressed 
by means of an accurate analysis, a simple estimate of the ways to minimize the energy of the system through the 
condensation of atoms at a relatively small number of the lattice sites suggests that patterns with higher values of 
the period may be favored by a still stronger contrast between the strengths of the DD and contact repulsion. An 
extension to multi-component condensates is interesting too. Further, Ref. 21], which recently reported a stable 2D 
checkerboard ground state in a related nonlocal model, suggests a possibility to search for stable density-modulation 



patterns in dipolar condensates trapped in a deep 2D OL, using the respective form of the DNLSE [41 1. 

On the other hand, while the analysis in this work is restricted to the DNLSE model based on the "frozen" lattice, 
the DD interactions between condensate droplets trapped at local minima of the OL potential may distort the droplet 
chain, making it necessary to add another degree of freedom — a shift of the droplets from local potential minima. In 
a combination with the equation for the droplets' mean-field wave function, this generalization may give rise to a new 
type of a dynamical system, blending the DNLSE with a model of the Frenkel-Kontorova type. Results obtained in 
the framework of such a coupled system will be reported elsewhere. Finally, it is interesting to extend the analysis to 
fermionic dipolar gases. In particular, it was very recently assumed that a Wigner crystal may emerge in a fermionic 
gas trapped in the OL [46j . 
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